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Abstract 

We apply the effective field theory approach to the three-nucleon system. In 
^ . particular, we consider S = 1/2 neutron-deuteron scattering and the triton. We 

^ \ show that in this channel a unique nonperturbative renormalization takes place 

which requires the introduction of a single three-body force at leading order. With 



one fitted parameter we find a good description of low-energy data. Invariance under 
the renormalization group explains some universal features of the three-nucleon 
C^ ' system — such as the Thomas and Efimov effects and the Phillips line — and the 

origin of SU{A) symmetry in nuclei. 
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1 Introduction 

Effective field theories (EFT's) are a powerful concept designed to explore a separation 
of scales in physical systems |]l| . For example if the momenta k of two particles are much 
smaller than the inverse range of their interaction 1/R, observables can be expanded in 
powers of kR. Following the early work of Weinberg and others [0], EFT's have become 
quite popular in nuclear physics f^, ^ . Their application, however, is complicated by the 
presence of shallow (quasi) bound states, which create a large scattering length a ^ R. 
In this fine-tuned case, a simple perturbative expansion in kR already breaks down at 
rather small momenta k ~ 1/a. In order to describe bound states with typical momenta 
k ~ 1/a, the range of the EFT has to be extended. A certain class of diagrams has to 
be resummed which generates a new expansion in kR where powers of ka are kept to all 
orders. For the two-nucleon system a power counting that incorporates this resummation 
has been found recently [^, §, 0]. A number of two-nucleon observables have been studied 
in the low-energy theory where pions have been integrated out Pj. Furthermore, pions 
can be included explicitly and treated perturbatively in this scheme P]. This approach 
has successfully been applied to A^A^ scattering and deuteron physics 0. 

The three-nucleon system is a natural test-ground for the understanding of nuclear 
forces that has been reached in the two-nucleon system. Furthermore, it shows some 
remarkable universal features. In the context of potential models, it has been found that 
different models of the two-nucleon interaction that are fitted to the same low-energy two- 
nucleon (NN) data predict different but correlated values of the triton binding energy 
-B3 and the S'-wave nucleon-deuteron (Nd) scattering length Og in the spin S = 1/2 
channel; all models fall on a line in the S3 x 03 plane, the Phillips line |l^ . Other 



universal features of three-body systems are the existence of a logarithmic spectrum of 
bound states that accumulates at zero energy as the two-particle scattering length 02 
increases (the Efimov effect fill), and the collapse of the deepest bound state when the 
range of the two-body interaction r2 goes to zero (the Thomas effect |[12|). These universal 
features have no obvious explanation in the framework of conventional potential models. 
(Alternative explanations of these features that are likely equivalent to ours can be found 
inRef. [0.) 

The critical issue in the extension of any theory from the two- to the three-body sector is 
the relative size of three-body forces. A naive dimensional analysis argument that does not 
incorporate the fine tuning in the two-body sector suggests that three-body forces appear 
only in higher orders in the kR expansion. We have shown that this is supported by Nd 



data in the S = 3/2 channel WM. However, for a system of three bosons a non-trivial 



renormalization enhances the size of three-body forces [0 . (For a different approach, see 



Ref. [|T6|.) Here we will show that the renormalization of the three-nucleon system in 
the S* = 1/2 channel is nonperturbative and requires a single three-body force at leading 
order. With such a force, the EFT allows to understand the above mentioned universal 
features in a unified way. Being SU{A) symmetric, this force also supports the view of 



spin-flavor as an approximate symmetry of nuclei [17]. Moreover, if its single parameter 



is fitted to 03 , we find a good description of the energy dependence of scattering and 



of triton binding. 

2 Preliminaries 

In order to avoid the difficulties due to the long-range Coulomb force, we concentrate 
here on the neutron-deuteron [nd) system. For simplicity, we also restrict ourselves to 
scattering below the deuteron breakup threshold where S*- waves are dominant. There 
is only one low-energy scale, \fJAB2 ~ 40 MeV, where B2 is the binding energy of the 
deuteron and M the nucleon mass. Since \fJAB2 is small compared to the pion mass, 
the pions can be integrated out, only nucleons remaining in the EFT as explicit degrees 
of freedom. The leading order of this pionless theory 0, |^, |^ is equivalent to the leading 
order in the KSW counting scheme [^ where pions are included perturbatively. As a 
consequence, the extension of our results to higher energies is well defined. 

There are two S'-wave channels for neutron-deuteron scattering, corresponding to total 
spin S = 3/2 and S = 1/2. For scattering in the S = 3/2 channel all spins are aligned and 
the two-nucleon interactions are only in the '^5'i partial wave. The two-body interaction 
is attractive but the Pauli principle forbids the three nucleons to be at the same point in 
space. As a consequence, this channel is insensitive to short-distance physics and precise 



predictions are obtained in a straightforward way |T^. There is also no three-body bound 
state in this channel. The S = 1/2 channel is more complicated. The two-nucleon 
interaction can take place either in the '^5'i or in the ^Sq partial waves. This leads to an 
attractive interaction which sustains a three-body bound state, the triton. The 5* = 1/2 
channel also shows a strong sensitivity to short-distance physics as the Pauli principle 
does not apply. We will see that the generic features of this channel are thus very similar 
to the system of three spinless bosons. There is in the latter a strong cutoff dependence 
even though all Feynman diagrams are finite, and renormalization requires a leading-order 



three-body force counterterm [15 



Let us start from the assumption that the three-body force is of natural size. The 
lowest-order effective Lagrangian is then given by 

£ = ^M ^-^o + ^ J AT - C* [N^T2acr2Ny ■ (AT^raaaaAr) (1) 

~C^ (N^a2TT2Ny ■ [N'^a2TT2N) + . . . , 

where the dots represent higher-order terms suppressed by derivatives and more nucleon 
fields, a (r) are Pauli matrices operating in spin (isospin) space, respectively. The 
contact terms proportional to Cq (Cq) correspond to two-nucleon interactions in the ^Si 
{^So) NN channels. Their renormalized values are related to the corresponding two-body 
scattering lengths a\ and 03 by Cg' = vrag /2M. Since no derivative interactions appear 
at this order, this Lagrangian generates only two-nucleon interactions of zero range. For 
practical purposes, it is convenient to rewrite this theory by introducing "dibaryon" fields 
with the quantum numbers of two nucleons fl^. In our case, we need two dibaryon 



fields: (i) a field T with spin (isospin) 1 (0) representing two nucleons interacting in 
the ^5*1 channel (the deuteron) and (ii) a field S with spin (isospin) (1) representing 
two nucleons interacting in the ^5*0 channel. Using a Gaussian path integration, it is 
straightforward to show that the Lagrangian (|ip is equivalent to 

C = Nnido + ^jN + ATT^-f+AsS^-S (2) 

-3T ^^t . N''T2aa2N + h.c.) - y (5^ ■ AT^aarrsAT + h.c.) +... 

At first it may look like the Lagrangian (^ contains more parameters than the original 
one, Eq. (|lD. However, the scales A^ and As are arbitrary and included in Eq. (|^) only 
to give the dibaryon fields the usual mass dimension of a heavy field. They can easily 
be removed by rescaling the dibaryon fields. All observables depend only on the ratios 
g^ g/Aj^^s- The introduction of these dibaryon fields makes explicit the formal similarity 



to the Amado model |]T9|. We stress, however, that the introduction of "quasi-particle" 
fields T and S carries here no dynamical assumptions 0. 

Since the theory is nonrelativistic, all particles propagate forward in time, the nucleon 
tadpoles vanish, and the propagator for the nucleon fields is 

Po — p /2M + le 

The dibaryon propagators are more complicated because of the coupling to two-nucleon 
states. The bare dibaryon propagator is simply a constant, i/As^r, but the full propagator 
gets dressed by nucleon loops to all orders as illustrated in Fig. |I]. The nucleon-loop 



Figure 1: Dressing of the bare dibaryon propagator. 

integral has a linear ultraviolet (UV) divergence which can be absorbed in g^ g/ At,s-, a 
finite piece determined by the unitarity cut, and subleading terms that have already been 
omitted in Eq. (^). Summing the resulting geometric series leads to 

-A5,T + -^^-Mp^ + 1^ _ ,e + ze 

where qt.s and At^s now denote the renormalized parameters. The A^A^ scattering am- 
plitude in the respective channel is obtained by attaching external nucleon lines to the 



dressed propagator. In the center-of-mass frame, the 5- wave amphtude in the ^Sq, ^5*1 
channels for energy E = k'^/M is 

^ ' M\ Mglr j 

and the renormahzed parameters gT,s and At,5 can be determined from 



(5) 



27rA 



S,T 



(6) 



Note that in this zero-range approximation, 



IMB2 = 1/4, (7) 

which holds only approximately, the discrepancy coming from range corrections. 



3 S = 1/2 nd Scattering and the Triton 

We now use the Lagrangian, Eq. (^, to describe nrf-scattering in the S = 1/2 channel and 
the triton. The two-nucleon interactions can take place both in the ^5*1 and ^Sq partial 
waves. There are two coupled amplitudes, a and b, as the triton can be built by adding 
a neutron to a proton and a neutron which are in either of the two partial waves. The 
amplitude a which has both an incoming and outgoing dibaryon field T gives the phase 
shifts for S = 1/2 nd scattering; a is coupled to the amplitude b which has an incoming 
dibaryon S and an outgoing dibaryon T. Although only a corresponds to elastic S = 1/2 
nd scattering, both amplitudes have the quantum numbers of the triton. The diagrams 
for the two amplitudes are obtained from the Lagrangian (|2|). However, the leading piece 
of all diagrams for a (6) is of order Mg^/Q"^ [MgTgs/Q"^), where Q stands for either the 
external momentum k or 1/02' . Therefore the diagrams have to be summed to all orders, 
which is conveniently done by solving the coupled integral equations shown in Fig. S. 



^^^ 



Figure 2: Coupled integral equations for S 
indicated by grey shaded (black) thick line. 



1/2 nd scattering. T (S) dibaryon is 



After the integration over the time component of the loop-momentum and the projec- 
tion onto the S*- waves has been carried out, we have 



-1 



^(^l/4 + v/3pV4-MEJ a{p,k) 

J2 



2 /-A q^ dq 
3p74 - ME - 1/a^ 



^(P' ^) + - / /,/ . K{p, q) [a(g, k) + 36(g, 
TT JO q^ — k'^ — le 



3K{p, k) + - f"^ -^^^—K{p, q) [3a(g, k) + h{q, k)] . 
IT Jo q"^ — k"^ — It 



2 i.2 HP'^) (9) 

2 /-A g2 ^^ 

Here k (p) denote the incoming (outgoing) momenta in the center-of-mass frame, ME = 
3k'^/4 — (1/02)^ is the total energy, and the kernel K{p, q) is given by 

T^( \ 1 , ( q^+pq-p^-ME \ 

K{p, q) = ^—ln{ — - , 10 

2pq \q'^ — pq — p^ — ME J 

The amplitude a(p, k) is normalized such that a(A;, k) = {k cot 6 — ik)~^ with 5 the elastic 
scattering phase shift. Furthermore, we have introduced a momentum cutoff A in the 
integral equations. Eqs. (||, |) have previously been derived using a different method pO 



In the limit A ^ 00 these equations do not have a unique solution because the phase of 
the asymptotic solution is undetermined ETI]. This is the same feature exhibited by the 



three-boson systems |]T5|. For a finite A this phase is fixed and the solution is unique. 
However, the equations with a cutoff display a strong cutoff dependence that does not 
appear in any order in perturbation theory. The amplitude a{p,k = const.) shows a 
strongly oscillating behavior. Varying the cutoff A slightly changes the asymptotic phase 
by a number of 0{1) and results in large changes of the amplitude at the on-shell point 
a{p = k). This cutoff dependence is not created by divergent Feynman diagrams. It is a 
nonperturbative effect and appears although all individual diagrams are superficially UV 
finite. 

In order to understand this strong A dependence, it is useful to note that the integral 
equations (§, |^) are SU{4) symmetric in the UV. Therefore it is sufficient to consider the 
SU{4:) limit {a\ = '^2 = '^2), since the cutoff dependence is a problem rooted in the UV 
behavior of the amplitudes. Furthermore, we note that the equations for a_|_ = [a + b] and 
a_ = [a — b] decouple in the S'f/(4) limit. The equations for a+ and a_ are 



-1 



^[l/a2 + ^j3pyA-MEJ a+{p,k) (li; 

-,2 



2K{p, k) + -f / f . 2K{p, q)a^{q, k) 
IT Jo q^^ — k^ — It 



^1 



\ (l/a2 + V3pV4 - MEJ a^{p,k) (12) 

2 /"A o^ dq 

= -K{p,k) / — -K{p,q)a_{q,k). 

IT Jo q^ — k^ — It 



The equation for a_|_ is exactly the same equation as in the case of spinless bosons while 
the equation for a_ is the same equation as in the S = 3/2 channel. The S = 3/2 equation 



is well behaved and its solution is very insensitive to the cutoff |Tj]. Consequently, the 
cutoff dependence in Eqs. (^, ^) stems solely from the equation for a+. As an example, 
the cutoff dependence of a+{p,k = 0) is shown by the solid, dashed, and dash-dotted 
curves in Fig. § for three different cutoffs, A = 1.0, 2.0, 3.0 x lO^a^^. 
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Figure 3: Cutoff dependence of a+{p,k = 0). Solid, dashed, and dash-dotted curves 
are ior H = and A = 1.0, 2.0, 3.0 x lO^a^^, respectively. Dotted, short-dash-dotted, 
and short-dashed curves show the effect of the three-body force for A = lO^a^^ and 
H = —6.0, —2.5, —1.8, respectively. 



Dependence on the cutoff reflects incorrect renormalization, and thus intolerable depen- 
dence on higher-energy modes. We need to ensure that only low-energy modes contribute 
explicitly, otherwise the power counting used to order two-body interactions — and in 
particular to produce Eq. (|) in lowest order regardless of the regularization procedure — 
would not hold. Now, since the equation for a+ is the same as in the boson case, we al- 
ready know the solution to this problem: a one-parameter three-body force counterterm 
H{A)/A^ that runs with the cutoff A W^. We simply replace 



K{p, k) -^ K{p, k) + 



g(A) 

A2 



(13) 



in Eq. (|TT]). The dotted, short-dash-dotted, and short-dashed curves in Fig. |^ show 
the effect of the three-body force H{A) on a+{p,k = 0) for A = lO^a^^ and H = 
—6.0, —2.5, —1.8. It is clearly seen that the variation of H for a constant A has the 



same effect on the amplitude as varying the cutoff. Consequently, we can compensate 
the changes in the asymptotic phase when A is varied by adjusting the three-body force 
term appropriately. (A more rigorous discussion of the renormalization procedure can be 
found in Ref. [T^.) 



We can obtain an approximate expression for the running of H{A) from invariance 
under the renormalization group. Requiring that the equation for a+ does not change its 
form when the high momentum modes are integrated out, we find 



H{A) 



sin(goln(A/A^) - arctg(l/go)) 
sin(so ln(A/A^) + arctg(l/so)) 



(14) 



where Sq ~ 1.0064 and A^ is a dimensionful parameter that determines the asymptotic 
phase of the off-shell amplitude [|l5l. The running of the three-body force H{A) according 
to Eq. (|1^) is shown by the solid line in Fig. ^. The dots are obtained by adjusting H{A) 
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Figure 4: Running of H{A) for A* = 0.9 fm 
numerical solution of Eq. (|TT|) (dots). 



(a) from Eq. (p^) (solid line), (b) from 



such that the low-energy solution of Eq. (|llD remains unchanged when A is varied. The 
observed agreement provides a numerical justification for our procedure. The three-body 
force is periodic with H{An) = H{A) for A„ = A exp(n7r/so) ~ A(22.7)". In particular, the 
bare three-body force vanishes for a discrete set of cutoffs. Note, however, that invariance 
under continuous changes in the cutoff does require a non-vanishing bare three-body force. 
An important point should be stressed here. We have first renormalized the two- 
body subamplitude (the dibaryon propagator) and then inserted the result in the three- 
body equation. The loop appearing in this equation was then regulated by a cutoff and 
renormalized by the introduction of a three-body force. This is equivalent to using separate 



cutoffs A (for the three-body equation) and A' (for the dibaryon propagator) and taking 
A' ^ oo first (changing the two-body parameters g^ ^/At^s as a function of A' to keep the 
two-body observables fixed) and then performing the renormalization of the three-body 
equation discussed above. A more standard approach would be to use the same cutoff for 
the loops appearing in both the dibaryon propagator and the three-body equation. The 
fact that the approach used in this paper is legitimate becomes evident by noting that 
the details of the renormalization are changed from one approach to the other, but not 
the physical results. To see that, let us consider the changes due to the use of A = A' 
throughout. The two-body subamplitude would then be changed by terms suppressed by 
powers of Q/A, where Q represents the typical momentum flowing through the two-body 
subamplitude. When inserted in the three-body equation, these terms are important only 
for loop momenta (g in Eqs. (§, ^) of the order of A. In this region of integration other 
arbitrary choices like, e.g., the precise shape of the cutoff used (sharp or smooth) are 
also important. All these differences occurring at the cutoff scale are absorbed in the 
three-body force, which would result in a different A dependence than given in Eq. (|14|). 
The physical results for small momenta (p -C A), however, are unchanged up to terms 
suppressed by powers oi p/A. Those terms depend on our particular choice of regulators 
and also break symmetries like Galilean invariance. There are many other terms of the 
same order in p/A that were already discarded, since we have performed an expansion in 
powers of p. The error introduced by this omission, however, is of higher order than we 
are working at. 

Formally, the three-body force term is obtained by adding 

Cs = -'-^^ {g^^N^if.^yif.m (15) 



+lgT9s[N^{f-ay{S-T)N + h.c. 
+glN^iS-TyiS-T)N), 



to the Lagrangian (^. Eq. (|TB|) represents a contact three-body force written in terms 
of dibaryon and nucleon fields. Via a Gaussian path integration it is equivalent to a true 
three-nucleon force, 

£3= - ^ME<^(Jl^(N-'^r,na,N)HN^<T,a,N)(N^T,a,a,N) (16) 



I 



A2 V4A 
^1_^T 9s 



't 



3 2 At 2A 



s 



{N''T2aka2N)\N^akTiN){N^a2TiT2N) + h.c. 



-A 



+ Th {N^(y2rkr2N)\NhuTiN){N^a2TiT^N) 

By performing a Fierz rearrangement, it can then be shown that the three terms in Eq. 
([T6| ) are equivalent. As a consequence, there is only one three-body force which is also 
SU{A) symmetric. Therefore the choice of the three-body force, Eq. (|1^), is by no means 

8 



arbitrary. In fact, it is the only S-wave three-body force with no derivatives that can 
be written down. This can be seen by remembering that three fermions (antifermions) 
in a completely symmetric spatial wave need to be in a completely antisymmetric spin- 
isospin state which is a 4 (4) of S'f/(4). Consequently, the three-body force transforms as 
4®4=1©15, but only the singlet is separately invariant under the spin and isospin 
subgroups. The two-body parameters g\l 1\t and Qs/ ^s appear in Eq. (|l^) because 
the effect of the three-body force cannot be separated from the effect of the two-body 
interaction. Both are linked by the renormalization procedure. Naive power counting 
would suggest that the three-nucleon force scales with l/{Mm\). The three-body force 
from Eqs. ([I5| , [T6|) , however, is enhanced by the renormalization group flow by two powers 
of a2. Using Eq. (|^), it is found to scale as a2/(Mm^) which makes it leading order. H{K) 
contains one new dimensionful parameter. A*, which must be calculated from QCD or 
determined from experiment. 

Recently Mehen, Stewart, and Wise p2[ remarked that the two-nucleon amplitude (H) 
is approximately SU{4) symmetric for p ^ l/a2- They also noticed that the only S-wave 
four-nucleon force that can be written down is SU{A) symmetric. Furthermore, there are 
no contact interactions with more than four nucleons without derivatives because of the 
Pauli principle. It is also reasonable to assume that the low-energy dynamics of other 
nuclei is dominated by S-wave interactions, as is the case for the deuteron and the triton. 
Therefore our S'?7(4) symmetric three-body force together with the findings of Mehen et 
al. 



22| gives an explanation for the approximate 5'[/(4) symmetry [^ in nuclei. 



Now we are in position to solve the full equations for the broken SU{A) case. Including 
the three-body force from above into Eqs. (^, P), we have 



^ (1/4 + pp^A - ME^ ' a{p, k) = K{p, k) + ^^ 
2 /-A q2 dq 



(17) 



+ 



TT Jo q^ — k^ — ie 



K{p, q)[a{q, k) + 36(g, k)] + ^^[a(g, k) + 6(g, k)] 



A2 



/3pV4 -ME- l/al 2H(A) 

2 ^ „ r-. — Kp, k) = 3Kip, k) + ^ ' 



p2 _ ^2 

g^ dq 
TT Jo q^ — k"^ — it 



2 /-A 
+ - 



A2 
K{p, g)[3a(g, k) + 6(g, k)] + ^^[a(g, k) + 6(g, k)] 



We need one three-body datum to fix the three-body force parameter A^,. We choose the 
experimental value for the 5* = 1/2 nd scattering length, 03 = (0.65 ± 0.04) fm |p3| , 
and find A* = 0.9 fm~^. (For the special cutoffs with vanishing H(A), we recover the 
results of Ref. IQ.) Although one three-body datum is needed as input, the EFT has 



not lost its predictive power. We can still predict (i) the energy dependence of S' = 1/2 
nd scattering and (ii) the binding energy of the triton. 

The resulting energy dependence oi S = 1/2 nd scattering for three different cutoffs is 
shown in Fig. ^. It is clearly seen that the introduction of the three-body force renders 
the low-energy amplitude cutoff independent. The scattering length is reproduced exactly 
because it was used to fix A*. The agreement for finite momentum is at least encouraging. 
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Figure 5: Energy dependence for S = 1/2 nd scattering with three different cutoffs 
A = 1.9, 6.0, 11.6 fm" (hidden under the sohd curve) for A* = 0.9 fm" . Dashed curve 
gives estimate of range corrections. Data are from the phase-shift analysis of van Oers 
and Seagrave ||25| (dots) and a measurement by Dilg et al. |2^ (triangle). 



Our experience from the S = 3/2 channel is that the range corrections improve the 
agreement considerably [jl4[. The dashed curve in Fig. ^ gives a crude estimate of these 
corrections. In our calculations we take the deuteron binding energy B2 from experiment 
and determine 05 from Eq. (^. The dashed curve in Fig. |^ is obtained by taking the 
experimental value of 03 as input and leaving all other parameters unchanged. From the 
estimated size of the range corrections, we anticipate an improved agreement once these 
corrections are included. (One should also keep in mind that the experimental phase-shift 
analysis |]25[ does not give any error estimate.) 

The triton binding energy is obtained from the solution of the homogeneous equations 
corresponding to Eqs. ([1^, ^8|) for E = — -B3. In Fig. | we show the bound state spectrum 
as a function of the cutoff A for a particular A^,. The shallowest bound state is the triton. 
Its binding energy is cutoff independent as long as A ^ 1/02- However, as A is increased 
new deeper bound states appear whenever H{A) goes through a pole. These new bound 
states appear with infinite binding energy directly at the pole. As the cutoff is further 
increased, their binding energy decreases and becomes cutoff independent as well. The 
poles of H[A) can be parametrized as 



An = /(A*a2) exp(n7r/so)a2 ^ . 



(19) 



One counts n bound states for A„_i < A < A„. However, only the states between 
threshold and A ~ m^ ~ 1/R are within the range of the EFT. By solving Eq. (p!9D for 
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Figure 6: Three-nucleon bound state spectrum for A* 
state corresponds to the triton. 



0.9 fm . The shallowest bound 



the number of bound states n, we then obtain 



#,,.!£,„ (I), 0(l) 



(20) 



and recover the well-known Efimov effect [^. In the limit a2 ^ oo, an infinite number 
of shallow three-body bound states accumulates at threshold. That these bound states 
are shallow follows from the fact that the binding energy is naturally given in units of 

oo. Furthermore, we also recover the Thomas effect 11121 . 



l/(Ma2) which vanishes as 02 

In a hypothetical world where A ~ 1/ R — * 00, the range of the EFT increases and deeper 
and deeper physical bound states appear. Consequently, there is an infinitely deep bound 
state for A — > 00. However, for m^r ~ 1/R as in the real world the deep bound states 
are outside the range of the EFT and their presence does not influence the physics of the 
shallow ones. 

Moreover, the variation of the parameter A^, gives a natural explanation for the Phillips 
line [|^. The three-nucleon system has been studied with many phenomenological two- 
nucleon potentials where particles can have momenta as high as 1 GeV. Different phase- 
equivalent two-body potential models have different off-shell amplitudes, which are equiv- 
alent to different three-body forces. In some cases explicit three-body potentials have also 
been added. These models can be viewed as particular examples of high-energy dynam- 
ics. At low energies, they must all be equivalent to the EFT with particular values for its 
parameters. What we have shown is that in leading order they can only differ in the value 
of the single three-body parameter A^.. Varying A* then generates a line. In Fig. ^ we 
show the Phillips line obtained in the EFT compared with results from various potential 
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model calculations [^ and the experimental values for B^ and a^ . Our Phillips line 
is slightly below the one from the potential models. We expect this discrepancy to be 
reduced once range corrections are taken into account. 
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Figure 7: S = 1/2 nd scattering length 03 as function of the triton binding energy B3 
(Phillips line): EFT to leading order (solid line); potential models (dots); and experiment 



cross . 



The dynamics of QCD chooses a particular value of A^,, which up to higher order 
corrections is A^, = 0.9 fm~^. For this value we find -B3 = 8.0 MeV for the triton binding 
energy, to be compared with the experimental result B^^^ = 8.48 MeV, which is known to 
very high precision. The theory without pions seems to work for triton physics. However, 
to draw definite conclusions one has to include the range corrections for both the energy 
dependence and the triton binding energy. 



The three-body force (|lg) involves three nucleons at the same point (within the resolu- 
tion of the EFT) in a relative S wave. It cannot contribute to the S = 3/2 channel where 
the three nucleons spin in the same direction. Indeed, we have found that the quartet can 
be well predicted with two-nucleon input alone [|I?|. We stress that the three-body force 
that arises at low energies is not the same that appears at higher energies. As momentum 
increases, part of the three-body contact operator (|T^ may be resolved into iterations 
of the two-body potential. The strong renormalization of this operator at low energies 
is not necessarily in contradiction with the potential-model-based phenomenology where 
three-body forces are found to be small. 
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4 Conclusions 

We have studied the three-nucleon system using EFT methods. While for nd scattering in 



the 5* = 3/2 channel precise predictions are obtained in a straightforward way [jT4| , the S = 
1/2 channel is more complicated. It displays a strong cutoff dependence even though all 
individual diagrams are UV finite. In this channel a nonperturbative renormalization takes 
place similar to the case of spinless bosons. This renormalization group argument requires 
an SU{4) symmetric three-body force which is enhanced to leading order. Together with 
the recent observations of Mehen et al. [^ this gives an explanation for the approximate 
SU{4) symmetry in nuclei. Furthermore, we find that the Phillips line is a consequence 
of variations in the new dimensionful parameter A,,, which is introduced by the three- 
body force. A* is not determined by low-energy two-nucleon data alone and has to be 
determined from a three-body datum. 

This is to be contrasted to other approaches to the problem that model the two-body 
interaction at distances of the order of the effective range as, for instance, the Amado 



model [|I9|. There a form factor is introduced in the interaction between two nucleons 
whose parameters are chosen to reproduce not only the two-body scattering length but also 
the effective range. Since the effective range expansion is not to be trusted at momenta 
of the order of the inverse effective range, this corresponds, in our language, to taking 
the three-body force equal to zero and picking a particular form for the regulator (the 
effective range term in the two-body subamplitude acts as a cutoff). 

The leading order of the EFT gives a quantitative description of the triton binding 
energy and the energy dependence for S = 1/2 nd scattering. In the appropriate limits 
for the two-body parameters 02 and i?, we also recover the well known Thomas and 
Efimov effects. The theory shows the potential for a realistic description of the triton once 
range corrections are included. Immediate applications of the EFT include polarization 
observables in nd scattering and triton properties such as its charge form factor. The 
incorporation of the long-range Coulomb force would widen the possible applications 
considerably as pd scattering and the physics of ^He becomes accessible. Work on these 
extensions is in progress. 

Finally, since in leading order the pionless theory is equivalent to a theory with explicit 
(perturbative) pions [H, the success in the three-nucleon system opens the possibility of 
applying the EFT method to a large class of systems with three or more nucleons. 
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